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Abstract 

No-scale supersymmetry or gaugino mediation augmented with large negative Higgs 
soft masses at the input scale provides a simple solution to the super symmetric flavor 
problem while giving rise to a neutralino LSP. However, to obtain a neutralino LSP it is 
often necessary to have tachyonic input Higgs soft masses that can give rise to charge- 
and-color-breaking (CCB) minima and unbounded-from-below (UFB) directions in the 
low energy theory. We investigate the vacuum structure in these theories to determine 
when such problematic features are present. When the standard electroweak vacuum 
is only metastable, we compute its lifetime under vacuum tunneling. We find that 
vacuum metastability leads to severe restrictions on the parameter space for larger 
tan (3 ~ 30, while for smaller tan [3 ~ 10, only minor restrictions are found. Along the 
way, we derive an exact bounce solution for tunneling through an inverted parabolic 
potential. 
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1 Introduction 

Super symmetry (SUSY) is a well-motivated extension of the Standard Model (SM) that 
provides an explanation for the stability of the hierarchy between the weak scale and the 
Planck scale [1]. However, experiment excludes SUSY from being an exact symmetry at 
low energies. If supersymmetry is softly broken, it is possible to push the superpartner 
masses up enough to make the theory consistent with experiment while still preserving the 
electroweak-gravity hierarchy. For this to work, the masses of the SM superpartners should 
not be significantly larger than the electroweak scale, putting them within reach of upcoming 
particle collider experiments such as the LHC 

A major drawback of this scenario of softly-broken low-energy supersymmetry is that 
the supersymmetry breaking operators generically introduce many new sources of quark and 
lepton flavor mixing, leading to flavor- changing neutral currents (FCNC) in conflict with 
experiment [2]. One way to cure this problem is to take No-Scale boundary conditions in 
which the soft scalar masses and couplings all vanish at a boundary input scale, with the 
gaugino soft masses being the only significant source of supersymmetry breaking in the visible 
sector at this scale [3]. Such boundary conditions can arise from gaugino mediation [4J, or 
conformal running effects [5], El [TJ El [9]. Low-scale scalar soft terms, which are necessary 
to obtain a viable phenomenology, are generated primarily from the gaugino masses in 
the course of renormalization group (RG) running from the boundary input scale down to 
near the electroweak scale. This dominant contribution is controlled by gauge interactions, 
implying that the regenerated scalar soft terms, and thus the squark and slepton masses, 
will be approximately flavor-universal at the low scale. 

A slight generalization of this pure No-Scale picture, which still leads to scalar soft terms 
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that are roughly flavor universal at the low scale (for the first and second generations), is 
to allow non-vanishing Higgs soft scalar masses at the boundary scale. In the context of 
gaugino mediation this can arise if the Higgs multiplets propagate in the bulk [5j, while for 
conformal running, it can emerge from hidden-sector interactions related to the origin of 
the fj, term [HI [U [10]. These Higgs Exempt No-Scale (HENS) models often have an impor- 
tant phenomenological advantage over pure No-Scale constructions: in large regions of the 
parameter space HENS models have a neutralino lightest superpartner particle (LSP) 
This differs from the situation in strict No-Scale models where typically the predominantly 
right-handed stau is the lightest SM superpartner particle \121. While such models can be 
viable if the true LSP is a gravitino [I3], a neutralino LSP with exactly conserved i?-parity 
tends to make for a better thermal dark matter candidate. 

HENS models can push the right-handed stau mass above that of the lightest neutralino 
by the influence of the non- vanishing Higgs soft masses on the RG running of the other soft 
parameters. The primary contribution to the slepton soft terms comes from the hypercharge 
D term Sy = Tr{Y m?) ~ ('^^f^ ~''^l/^)/2 that appears in the RG equations for all the 
MSSM scalar soft masses [lll[T5l[T2]. This contribution increases the right-handed slepton 
soft masses in the course of running down provided Sy is negative. In the portion of the 
HENS model parameter space that is also consistent with electroweak symmetry breaking, 
the Higgs scalar soft masses must typically be very tachyonic at the boundary scale to obtain 
a neutralino LSP, and the uncolored sfermions still end up being relatively light [llj. 

While these HENS models are compelling, the tachyonic Higgs soft masses at the bound- 
ary scale lead to concerningly large tachyonic low-scale Higgs soft masses. These, combined 
with the small (but positive) slepton soft masses, suggest that the true vacuum of the theory 
might be a charge-and-color-breaking (CCB) minimum, or there may exist unbounded-from- 
below (UFB) directions that are only stabilized far out in field space by higher-dimensional 
operators [161 [T71 [HI [191 [201 |2l| . The presence of such features need not exclude these regions 
of the parameter space provided our SM electroweak vacuum is metastable and long-lived 
relative to the age of the universe. 

In the present work we study the vacuum structure of HENS models, and compute the 
lifetime of the SM vacuum state when it turns out to be metastable. We concentrate on the 
stability of the SM vacuum at zero temperature with respect to vacuum tunneling. Thermal 
effects in the early universe can also potentially induce thermal transitions between different 
vacua, and modify the shape of the potential itself. However, these thermal effects tend 
to stabilize the origin of the field space to which the SM vacuum is connected, favoring 
this vacuum over others that lie further out in the field space [22|. Thermal effects are 
especially effective in delaying the formation of vacua that break color due to the large 
thermal corrections from the strong gauge and top quark Yukawa couplings [22l [201 [23] . 
Scalar field evolution during and after inflation may also populate non-standard vacua, 
although the precise result depends on the details of inflation and lies beyond the scope of 
this paper [2H [21]. By focusing solely on the T = constraints we obtain conservative and 
unambiguous bounds on the HENS parameter space that do not depend on the cosmological 
history. 

Motivated by gauge coupling unification, throughout this work we take the input scale 
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to be Mc = 10^^ GeV and we assume a universal gaugino mass mi/2 at this scale. The set 
of independent free parameters after imposing consistent electroweak symmetry breaking is 
therefor^ 

mi/2, tan/3, m^^jM,), m]j^{M,). (1) 

Investigations of No-Scale models with non-universal gaugino masses can be found in Refs. [25| 
[26] . Similar input soft terms have also been studied in the context of Refs. [271 128] 

The outline of this paper is as follows. In Section [2] we investigate the vacuum structure 
of HENS models. We discuss the general features of vacuum tunneling and outline the 
relevant criteria we use to judge which of the possible CCB and UFB vacua are the most 
dangerous in Section [31 In Section [H we map out the portions of parameter space that are 
phenomenologically consistent and give rise to a SM vacuum state that is sufficiently long- 
lived. Section [51 is reserved for our conclusions. In Appendix [Al we outline our numerical 
technique for estimating the vacuum tunneling rate. We present in Appendix [Bl an exact 
tunneling solution for an inverted-parabola barrier. 



2 UFB Directions and CCB Minima in the HENS Model 

With the HENS model input parameters of Eq. ([T]), relatively small slepton soft masses 
as well as large tachyonic Higgs soft masses obtain near the electroweak scale in much of 
the phenomenologically allowed parameter space [11]. This spectrum of soft parameters 
frequently implies the existence of UFB directions and CCB minima that are much deeper 
than the standard electroweak vacuum. Indeed, we find that nearly the entire allowed 
parameter space in the HENS model (before imposing vacuum stability constraints) has 
at least one UFB direction. We investigate the existence and nature of such potentially 
dangerous vacuum features in the present section. 



Given the soft breaking spectrum that arises in the HENS scenario, the results of Ref. [18] 
suggest that the most dangerous vacuum feature will be a sleptonic UFB-3 direction. This 
direction has Hd = 0, with fi, fn, i^ii^^, and all non-zero. Turning on expectation 
values for these fields, there exists a D- and F-fiat direction that is only lifted by quadratic 
supersymmetry breaking operators. To obtain F-fiatness, only the F-term of Hd must be 
cancelled. This can be arranged by taking 



Vr 



(2) 



with the relative phase of tl and r/j chosen appropriately. D-fiatness is then obtained by 
setting 

K«l^ = -(^) + |f|^ + l^^-T. (3) 



-'^The sign of the /i term is also free, but we fix it to be positive. This shghtly increases the model 
prediction for (g — 2)^, as favored by the experimentally measured value of this quantity 
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and represents the lowest-energy F-flat field configuration provided 



m > 



+ 



4m|. 
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'2yr 



(4) 



The scalar potential along this direction in field space then becomes [18 



Vj 



UFB-3 



2m\. 

g^ + g 



(5) 



When {rn\i^ + Jti^.) is negative, the potential becomes unbounded in the limit — > oo. 
It will ultimately be stabilized by loop corrections or higher-dimensional operators (that we 
have not included in Eq. ([5])) at a location that is very deep and far out in field space relative 
to the electroweak vacuum. 

This sleptonic UFB-3 direction is particularly dangerous in the HENS models on account 
of the large and negative values of m|^^ and the smaller values of m|. and m|;, that emerge 
in the low-energy spectrum. These properties imply that the barrier against tunneling from 
the electroweak vacuum near the origin out to the deeper UFB-3 direction, arising from 
the linear term in in Eq. ([5]), will not be especially large. The barrier will be further 
weakened by larger values of tan/3 which enhance the coupling = irir/v cos j3. Other 
similar UFB-3 directions may be present in the theory, but they will generally have larger 
barriers due to the larger values of the squark soft masses or the smaller values of the first- 
and second-generation lepton Yukawa couplings. 

When \H^\ does not satisfy the bound given in Eq. (jl]), the lowest-energy F-fiat direction 



in the potential has ji^Li^al = 0, and is given by [18 
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This potential is no longer D-fiat, and is stabilized at a finite value of \H^\. If this point 
occurs with less than the bound of Eq. (jlj), it is a constrained local CCB minimum. 
On the other hand, when this constrained local extremum has \H^\ larger than the bound 
of Eq. (jl]), it represents a saddle point that is unstable under flowing to a non-zero value of 

In practice, we find that for smaller values of \H^\ the potential can be reduced further by 
relaxing the F-fiatness constraint of Eq. ([2]). The effect of dropping the F-fiatness constraint 
is that the minimal potential for a given value of \H^\ is deformed slightly away from the 
precise UFB-3 form of Eqs. f l5ll6l) . but that the same qualitative features remain. In particular, 
there usually remains a local extremum at ^ and l^'L^.^gl = 0. We shall designate this 
local extremum as CCB-4. If this extremum occurs at smaller values of \H'^\, on the order of 
the electroweak scale, it can be a local CCB minimum. When the CCB-4 extremum occurs 
with a value of much larger than the electroweak scale, it is generally a saddle point 
that flows in the v^.^^ direction to a genuine UFB-3 direction of the form given in Eq. ([5]). 
Even when it is only a saddle, the CCB-4 point plays an important role in determining the 
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tunneling rate from the electroweak minimum to the UFB-3 direction, as we will discuss 
below. 

A stable CCB-4 minimum with stops can also arise for smaller values of tan (3 and a 
correspondingly larger yt Yukawa coupling. In this case we have = while I^lI = \tn\ = \t\ 
and are all non-zero. The relevant potential is 

Vd = m^Jif,T + K + Ol^T + \f^H,- y,P\^ + {\H,\^ + . (7) 

This potential is generally stable against excursions in the \d\ = \dL-^.^\ = \dR-^^\ direction 
on account of the larger squark soft masses that arise in the HENS model. With m^^ < at 
the low scale, this potential often has a CCB-4 minimum with both \H^\ and \i\ non-zero. 
However, as we show below, the barrier to tunneling to this minimum from the standard 
electroweak minimum to this CCB-4 minimum is usually safely large, again on account of 
the larger values of the squark soft masses as well as the less negative values of m'jf^. For 
similar reasons, we expect that the rate for tunneling to the other potential CCB minima 
discussed in Ref. [TS] will typically be less constraining than the rate to tunnel to a stau 
UFB-3 direction. 



3 Computing the Vacuum tunneling Rate 

The existence of vacua deeper than the standard electroweak minimum in HENS models 
implies there is a danger of tunneling into one of these phenomenologically unacceptable 
states. At the very least, the lifetime for this tunneling must be greater than the age of the 
universe. The vacuum-to-vacuum transition rate associated with tunneling can be calculated 
using path integral methods [29l [30j . In the semiclassical approximation, the lifetime of the 
vacuum is found to be 

J- = r/V = Ae-^^^^^ (8) 

TV 

where A is a dimension-four prefactor to be discussed below, (pi denotes the bounce solution 
for the i-th field, and Sb is the Euclidean action, 

SS] = j d^xs(|V0,|' + f/(0i,...,0,)) =T[^^] + V[U. (9) 

The bounce solution for the fields is the extremum of the Euclidean action that is 0(4)- 
symmetric and obeys the the following boundary conditions: 

^ ^ ' =0 (10) 



dp 

hm Up) = (11) 



p— >oo 



where p = (x^ + t'^Y^'^ is the Euclidean distance. These boundary conditions correspond to 
a field that begins in a metastable vacuum and tunnels through a barrier separating it from 
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a deeper vacuum state, emerging with zero kinetic energy. We focus on 0(4) symmetric 
solutions because these are expected to have the least action, and therefore dominate the 
tunneling probability [31]. 

Even with the simplification of an 0(4) symmetry, the equation of motion for the bounce 
cannot in general be solved analytically. We compute the bounce solution numerically using 
the improved action method ^32j . The details of this method are outlined in Appendix lAl 
Even more difficult to compute is the non-exponential pre-factor A in Eq. ([8]) [30j. On general 
grounds, we expect it to be on the order of the mass scale setting the size of the potential 
barrier. In low-energy supersymmetry, a good estimate for this number is (100 GeV)^, which 
we take to be the case throughout the rest of this article. The precise value of this pre-factor 
is unlikely to affect our qualitative conclusions as the multiplicative uncertainty in its value 
is much more slowly-varying than the exponentiated large values of the action that lead to 
acceptable lifetimes. Our choice for the pre-factor is also conservative, in that choosing a 
larger number here would only exclude more points. With this pre-factor, it is found that 
the lifetime of the SM vacuum will be greater than the age of the universe, l/(r/V^) > t^, 
provided Sb[4>] > 400. 

The value of the bounce action for tunneling between a pair of local vacua depends on 
the relative depth of the minima, the number of distinct minima, the height of the barriers 
between them, and the relative size of the field values within them. In general, the bounce 
solution represents a configuration of fields that simultaneously minimizes these opposing 
contributions, with the kinetic term favoring slowly varying fields, and the potential term 
preferring to reach the deepest minimum as quickly as possible. Some intuition for this 
balance can be obtained by examining special cases in which the bounce solution can 
be derived analytically. One such example is a potential consisting of piecewise linear 
segments [M]. Motivated by the UFB-3 potential in Eq. ([5]), we present another exact 
tunneling solution in Appendix [B] for an inverted parabolic potentialo 

As expected, our analytic tunneling solution shows that the bounce action increases with 
the size of the barrier. This solution also implies that the depth of the minimum to which 
one is tunneling to ceases to matter once it becomes very deep. Thus, we can safely compute 
the rate to tunnel into a UFB-3 direction without knowing where it is ultimately stabilized 
provided the corresponding minimum is very deep relative to the height of the barrier. On 
the other hand, the relative depth of the minima is relevant to the tunneling rate when 
they are nearly degenerate, as can be seen from the analytic thin-wall approximate solution 
that can be safely applied in this case |29]. One important feature not captured by our 
one-dimensional analytic tunneling solution is that the bounce action tends to also increase 
when there are more independent fields involved in the tunneling process, as each one of them 
contributes non-negatively to the kinetic portion of the bounce action. For these reasons, 
numerical solutions are used to get the details correct. 

^ While this potential has the same form as the UFB-3 potential given in Eq. (O, we cannot apply our 
exact solution to tunneling through the UFB-3 barrier since the relevant kinetic terms for this case are 
non-canonical. 
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Figure 1: The bounce action for tunneling to a stau UFB-3 direction or a CCB-4 minimum as 
a function of the HENS model high scale input parameters m|^^ (Mc) and (Mc) for mi/2 = 
500 GeV and tan (3 = 30. All points shown are consistent with collider phenomenology. The 
points enclosed below by the dash-dot line have a neutralino LSP. The solid line separates 
the region with a CCB-4 minimum or saddle point (left) from that which only has a UFB-3 
(right). 5* > 400 is cosmologically safe. 



4 Vacuum Stability Bounds on the HENS Parameters 

The possibility of tunneling from the SM vacuum to a phenomenologically unacceptable 
vacuum places strong constraints on the parameter space of the HENS model. To investigate 
these constraints, we have calculated the bounce action for tunneling to a CCB-4 minimum 
or a UFB-3 direction for a series of representative HENS parameter sets. Our strategy is to 
fix mi/2 and tan/3, and scan over the input values of m^j^^M^ and m\^{Mc) that lead to an 
acceptable low energy spectrum. We focus on the values mi/2 = 500, 1000 GeV for tan /5 = 30 
and mi/2 = 300, 500 GeV for tan /5 = 10. In each of these scans we have consider points that 
meet the phenomenological constraints laid out in Ref. ^JJ, including the current LEP lower 
bounds on the lightest Higgs boson. We also keep the points with a charged slepton as 
the lightest MSSM superpartner found in Ref. [TT], noting that they would require a more 
complicated cosmology and likely R-parity violation or a gravitino LSP [13]. By keeping 
such points, our results are also applicable to minimal gaugino mediation subject to our 
assumptions about the gaugino masses and the input (compactification) scale. 

We first consider the cases with tan (3 = 30, where the tau Yukawa coupling is enhanced. 
This has the effect of opening up a CCB-4 minimum or CCB-4-like saddle point that flows 
to a UFB-3 direction. In Fig. [T] we show ranges of the bounce action for tunneling out of the 
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Figure 2: The bounce action for tunneling to a stau UFB-3 direction or a CCB-4 minimum 
as a function of the high-scale HENS model input parameters m|^^(Mc) and m'jj (Mc) for 
mi/2 = 1000 GeV, tan/? = 30 and sgn{fi) = 1. All points shown are consistent with collider 
phenomenology. The points enclosed below by the dash-dot line have a neutralino LSP. The 
solid line separates the region with a CCB-4 minimum or saddle point (left) from that which 
only has a UFB-3 direction. S > 400 is cosmologically safe. 



SM vacuum to a CCB-4 minimum or a UFB-3 direction for mi/2 = 500 GeV as a function 
of the input values of m'\j^{M^ and m^^{Mc). All points to the left of the solid line in this 
figure have either a CCB-4 minimum or saddle point. The points to the right have either a 
stable SM vacuum or a UFB-3 direction with no CCB-4 stationary point. To the left and the 
right of the solid line we show the bounce action for tunneling along the direction leading to 
the CCB-4 stationary point or to the UFB-3 direction, whichever is smallerjl We emphasize 
the CCB-4 stationary point here, even when it is only a saddle point, because we find that 
the most dangerous lowest-action tunneling path from the electroweak minimum is typically 
one that passes near this point with f'i.^g = 0. The dot-dashed line in Fig. [1] indicates the 
upper border of the portion of parameter space in which the LSP is the lightest neutralino. 
(See Ref. [11] for more details.). In Fig. [2] we show the same quantities as in Fig. [1] for 
mi/2 = 1000 GeV. 

Both Fig. [1] and Fig. [2] show that for tan f3 = 30 the lifetime for tunneling to a stau 
UFB-3 direction or CCB-4 minimum is shorter than the age of the universe, corresponding 
to Sb < 400, over a large portion of the otherwise acceptable parameter space. As expected, 
the newly disallowed regions are those with large and negative values of and m|^^ at 

■^Except for those points very near the sohd hne, the CCB-4 extremum is a saddle point flowing to the 
UFB-3 direction. 



9 




Figure 3: The bounce action for tunneling to a stau UFB-3 direction as a function of the 
HENS model parameters m|^^ and Tn'jj^. The other HENS parameters have been fixed to 
be mi/2 = 500 GeV, tan/5 = 10, and sgn{fi) = 1. All points shown are consistent with 
collider phenomenology. The points between the two dash-dotted lines have a neutralino 
LSP. 5* > 400 is cosmologically safe. 



the input scale M^. Of these two soft masses, only fnjj^ appears in the potential relevant 
for the UFB-3 direction or the CCB-4 minimum, and it has the stronger effect. Indeed, 
the isocontours of Sb coincide roughly with lines of constant fnjj^ for smaller values of fn'jj^. 
For smaller values of fn'jj^ the SM vacuum becomes sufficiently long-lived to describe our 
universe. The bounce action for these points is much larger simply because the effective 
width of the barrier is larger. Note that this stable region includes the minimal gaugino 
mediation point, fn'jj^ = Tn\j^ = at Mc- 

Comparing the plots for mi/2 = 500 GeV and mi/2 = 1000 GeV, we see that larger mi/2 
tends to yield a slightly more stable SM vacuum. This arises simply because increasing the 
input gaugino masses also increases the low-scale slepton soft masses through RG running. 
On the other hand, larger gaugino masses also permit more negative values of m|^^ , so there 
remain significant parameter regions in which the tunneling rate is too fast. From Fig. [H 
where mi/2 = 500 GeV, we see that nearly the entire region in which the lightest superpartner 
is a neutralino is ruled out by our vacuum stability considerations. For mi/2 = 1000 GeV, 
there is a small region in which the lightest superpartner is a neutralino and the electroweak 
vacuum is sufficiently long-lived. 

We turn next to the cases with tan/3 = 10, for which the tau Yukawa coupling is less 
enhanced. In Figs. |3] and H] we show ranges of the bounce action for tunneling from the SM 
vacuum to the UFB-3 direction for tan/3 = 10 and M1/2 = 500, 300 GeV as a function of 
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Figure 4: The bounce action for tunneling to a stau UFB-3 direction as a function of the 
HENS model parameters fn'jj^ and fn'jj^. The other HENS parameters have been fixed to be 
mi/2 = 300 GeV, tan/3 = 10, and sgn{fi) = 1. All points shown are consistent with collider 
phenomenology. The points between the two dash-dot lines have a neutralino LSP. S > 400 
is cosmologically safe. 



the input values of Tnjj^ and at M^. No CCB-4 local extremum is found for any of the 
points scanned over since the deviation from F-flatness is a sizeable effect in this case. The 
dot-dashed lines in these plots enclose the portion of parameter space in which the LSP is 
the lightest neutralino. (See Ref. [H] for more details.) From these plots we see that the 
lifetime of the standard electroweak vacuum against tunneling to a UFB-3 direction is safely 
large for both mi/2 = 500, 300 GeV. This is the result of the smaller value of the tau Yukawa 
coupling i/t-, which gives rise to a larger barrier against tunneling, Eq. ([5]). As before, with 
all else equal Sf, increases with mi/2 and so higher values of mi/2 remain safe. 

At smaller values of tan/? the top quark Yukawa coupling grows larger, and we should 
check that the tunneling rate to CCB-4 minima associated with the stops is adequately small. 
As was previously discussed, this direction will only occur if is large and negative. In 
Fig. [5] we show contours of the bounce action for tunneling to the stop CCB-4 minimum, as 
well as the regions in which the SM minimum is the true minimum. Only a very few points 
at the largest and most negative values of are excluded, while in the great majority 
of the parameter space the SM vacuum is the true minimum. We find a similar result for 
mi/2 = 300 GeV. 
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Figure 5: The bounce action for tunneling to a stop CCB-4 minimum as a function of the 
HENS model parameters m|^^ and 'm'jj^. The other HENS parameters are mi/2 = 500 GeV, 
tan/3 = 10 and sgn{fi) = 1. All these points are consistent with collider phenomenology. 
The points between the two dash-dot lines have a neutralino LSP. 5* > 400 is cosmologically 
safe. 



5 Conclusion 

In the present work we have examined the constraints placed on the HENS model from 
vacuum stability. Due to the large tachyonic Higgs soft masses that can emerge in this model, 
there often arise local vacuum states deeper than the standard electroweak minimum. Many 
points that are consistent with collider phenomenological constraints (described in Ref. [llj ) 
are ruled out because they lead to an overly short-lived SM vacuum. 

The most dangerous vacuum feature is a UFB-3 direction involving the stau fields. We 
have also found a new CCB-4 saddle point that facilitates tunneling to the UFB-3 direction. 
As a result, vacuum tunneling rates tend to be too fast for larger values of tan/5. At lower 
values of tan /3, tunneling to a CCB-4 direction involving stop fields rules out a very small 
portion of the parameter space, with the rest of the parameter space being safely long-lived. 

We conclude that the HENS models with a neutralino LSP and larger values of tan f3 = 30 
are mostly ruled out subject to our assumptions about the input scale and gaugino mass 
universality. On the other hand, minimal gaugino mediation and the HENS models without 
a neutralino LSP may still have a sufficiently long-lived electroweak vacuum state at larger 
values of tan/?. For lower values of tan/?, such as tan/5 = 10, the constraints from vacuum 
tunneling are much weaker, and most of the parameter space that is consistent with collider 
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lower-energy bounds remains viable. 
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A Appendix: The Improved Action Method 

The bounce action is a stationary point of the Euclidean action given in Eq. subject 
to the boundary conditions of Eqs. ( fTOfTTl) . The corresponding equations of motion for the 
0(4) symmetric solution are 

ap^ p dp ocpi 

where i runs over the independent fields. These equations are a set of non-linear coupled 
differential equations with an a priori unknown starting point. These conditions together 
make it a very difficult problem to solve and require numerical techniques [TTl [321 133] . 

The technique we use in the present work is called the improved action method |32j . 
In this method, additional terms are added to the action that are identically zero for the 
bounce solution. The advantage of adding these terms is that they make the bounce solution 
a minimum of this modified action and not just an extremum. The term that does this is 
found by making the change of variable x ^ ax in Eq. iQ. Because the bounce is the 
extremum of the action, the first derivative of the scaled action with respect to a will be 
zero for a = 1. This gives the following condition: 

{a^T[^,] + 2a^f/[0,,]) = (13) 

This relation illustrates that the potential term must be negative. The kinetic term cannot 
be negative because it is the integral of a sum of squares. Since the potential term scales as 
and the kinetic term scales as cr^, Eq. ( JT3l) defines a maximum. Thus, we have determined 
the maximal direction of the saddle point. By adding to the action the absolute value of this 
quantity to a positive power, the saddle point of the action can be turned into a minimum. 
In this case the improved action is 

^[0,,] = r[0i] + U[^i] + A |T[0,] + 2t/[0,]|" (14) 

where A and n are positive constants. 
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To solve for the bounce with this improved action, we take an initial profile for the vevs 
0j with the kinetic and potential terms 

m=l 

L-1 

U[^,] = 27r'A'Y,iPm+l- Pm)piU{^T,-An)- (16) 
m=l 

A is a parameter determined by Eq. ( fT3l) . Inspired by the thin- wall approximation, we take 
the following initial guess for the bounce solution 

0(p) = atanh(6(p — po)) + c. (17) 

The coefficients a and c can be solved for by applying the boundary conditions 0(0) = 0e 
and 0(00) = (pf, where (pj are the field values in the SM minimum and 0e are the values in 
the vacuum to which the tunneling connects. Since 0e is a priori unknown, this leaves (pe, b, 
and po as free parameters. These parameters are determined by first substituting the field 
profile in Eq. (fT7|) into the modified action for each 0j. 

The 0e will not in general be the field configuration of the minima, but rather some field 
points inside the well. In the case of a UFB-3 direction there is no minimum, and 0e will be 
some point on the runaway downslope with a potential energy less than that at SM minimum. 
The exact value of 0e as well as b and po are determined using a minimization routine that 
finds the coefficients that minimize the modified action. In Ref. [20] the authors used the 
thick-wall approximation as an initial guess. In our case, the guess given in Eq. (fT3!) works 
better numerically because it is adaptable to both thick- and thin-wall potential profiles. 
Once the initial profile, coefficients and all, is determined, we randomly vary each lattice 
site. The variations are stopped when further iterations do not reduce the modified action. 
To ensure that we arrive at the bounce solution, we choose a value of A that ensures 0.999 < 
(— T[0j]/2/\^[0j])^/^ < 1.001. The smallest value of A able to maintain this condition and 
used to optimize the code was close to 0.5 with n = 1. 



E 



I m+l 



Lm\2 



(15) 



B Appendix: Tunneling Through an Inverted Parabola 

We exhibit here an exact bounce solution for a single field tunneling through an inverted 
parabolic potential. To the best of our knowledge, this simple solution has not been presented 
elsewhere in the literature, although a related solution for a linearized potential was obtained 
in the pioneering work of Ref. 



Consider the run-away potential 

V^(0)=a|0|-6|0|^ (18) 

where is a real scalar field. The origin is a metastable minimum for a, b > 0. Focussing 
on > 0, the equation of motion for an 0(4)-symmetric bounce is 

_| + __2 + (2*^_„) = 0, (19) 
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where p is the Euchdean distance. The boundary conditions for the bounce are 



dp 



= 0) 
ip = R) 









(20) 
(21) 
(22) 



For potentials in which the false vacuum is flat, R oo, as we assumed in Eq. (fTTj) . Here, 
the potential has a singular first derivative at the false vacuum at the origin. As a result, 
the bounce reaches the origin at finite R and remains there asymptotically [31]. We have 
checked that deforming the potential into a concave parabola very near the origin to resolve 
this singularity leads to i? — > oo, but otherwise has only a small (and smooth) effect on the 
solution presented here away from the origin. 

Eq. (fT9|) is a linear second-order differential equation that can be solved analytically. The 
general solution before imposing boundary conditions is 



2,0 X 

where Ji and Yi are Bessel functions of the first kind, 

P 

X = 



2b 



(23) 



(24) 



and Ci and C2 are constant coefficients. The values of Ci and C2, as well as R, are fixed by 
the boundary conditions, Eqs. (12011221) . The first of these, Eq. (120|) . gives 







(-[CiJiix)+C2Yi{x)] 

ax \x 



x=0 



— lim 



X 



[cMx) + C2Y2{x)] 



(25) 



where we have made use of the recursion properties of Bessel functions [35j . For this limit 
to be non-singular, we must have C2 = 0. On the other hand, J2{x)/x vanishes at the origin 
so no constraint is imposed on Ci by this condition. Applying Eq. ( l22l) . we find 







d_ 

dx 



-Ji{x) 

X 



(26) 



-Ci J2(Xh 



where we have defined x+ = R/y2h. From this, we conclude that x+ must be a zero of J2. 
In particular, for the minimal action bounce solution, it should be the first non-trivial zero: 



x+ ^ 5.1356 . . . 



(27) 
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Finally, let us apply the condition of Eq. (12T]) to fix ci, 



ci = - 



2b 



_Ji{x+) 



{21 



We can use this exact solution to compute several interesting quantities related to the 
bounce. The bounce action is found to be 



271^ / dpp^ 
Jo 

32''+¥ 
(215)^. 



(29) 



The value of the field (p at the escape point where it emerges from tunneling is given by 



m 



a a 
26 ~ 26 



Ji(x . 



x^O X 



Aa 



(30) 



For comparison, the maximum of the potential occurs at = a/26 and the width of the 
potential barrier is A0 = a/6. Thus, the escape point is well beyond the peak, but only by 
a factor of a few. 

The exact solution presented here can also be extended to potentials that are piecewise 
segments of parabolas, in analogy to the solution for piecewise-linear potentials presented 
in Ref. [Mj. (When the local potential is a concave parabola, the general local solution 
consists of modified Bessel functions of the first kind.) This offers the possibility of obtaining 
closed-form expressions for bounce solutions to a wide variety of potentials by approximating 
them with segments of parabolas. In practice, we find that the matching conditions between 
adjacent parabolic segments often lead to complicated implicit transcendental equations. The 
generalization to thermal transitions over piecewise-parabolic barriers is also straightforward. 
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